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Abstract 

It has been known for many years that there is a CP-conserving component for 
the decay mode Ki 'K^e^e~ and that its magnitude can be obtained from a mea- 
surement of the amplitudes in the Kl tt^^^ decay mode. We point out that the 
usual description of the latter in terms of a single parameter, ay, is not sufficient to 
extract the former in a model independent manner. We further show that there ex- 
ist known physics contributions to t:^^^ that cannot be described in terms of 
the single parameter Uy. We conclude that a model independent analysis requires the 
experimental extraction of three parameters. 




1 Introduction 



The mode '7r°77 has been the subject of intense study both as a test of chiral pertur- 

bation theory and as the source of a CP-conserving amphtude for —>■ -K^e^e" 
It has been known since the first experimental results appeared [^] that lowest order {p'^) 
chiral perturbation theory is not sufficient to explain simultaneously the observed rate and 
spectrum. For some time now, it has become standard to use a theoretical description which 
incorporates certain non-analytic terms at next to leading order (p^) |]TD|, |TT[, as well as one 



parameter, ay |TT]]. This parameter arises in vector meson dominance models for this decay 



[|T2|, but it does not parameterize the most general analytic amplitude at next to leading 



order in chiral perturbation theory, p^. Instead, at order p^ the amplitude is described by 



three independent parameters: a^^, and [3 in the notation of Cohen et al. |Tl[] 



Nevertheless, the parameterization of the amplitudes for 7r°77 in terms of 

Qy alone has been retained in the literature. In this paper we wish to point out that 
this is insufficient if one wants to extract a model independent bound on the CP-conserving 
component of n^e^e^ from experiment. This is something which should be considered 

by the forthcoming experimental analyses of the mode 7r°77 by the KTeV and NA48 

collaborations. Within the framework of chiral perturbation theory the new data should be 
analyzed in terms of a^^, 1^- The issue of whether vector mesons dominate this decay 

mode is an experimental question, and should not be an input in the analysis of data. As a 
further motivation for the more general fit, we show in this paper that there exists known 
physics, the /2(1270), which affects the Kl 7r°77 amplitude at a level comparable to that 
of vector mesons, and which cannot be parameterized in terms of the single constant ay. It 
should be no surprise that the /2(1270) can play an important role in this decay mode, given 



its prominence in the reaction 77 tt^'tt' 



Of particular importance is the determination of the CP-conserving contribution to 
Kl —>■ 7r°e^e~. This contribution is completely dominated by one of the two amplitudes 
present in Kl —>■ 71^77. In fact, it is of phenomeno logical relevance only when it arises from 
the amplitude in which the two photons are in a relative D-wave For this reason an 
accurate determination of both amplitudes is crucial. The model independent analysis we 
advocate here permits the extraction of the necessary information directly from the data, 
whereas the usual analysis in terms of ay forces correlations between the two Kl —* 7r°77 
amplitudes which may or may not be present in the data. 



In this section we review the parameterization of the Kl ir^^yy amplitude with terms of 
order up to in chiral perturbation theory and we compare fits to the KTeV data from 
1999 in terms of ay and in the general parameterization. 

The most general form of the K TC'-f'-f amplitude contains four independent invariant 
amplitudes A, B, C and D and has been described in the literature before ||^. For the case 
of Kl —>■ TT^jj, and in the limit of CP conservation, only two of these amplitudes come into 
play: 
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where Gg = 9.1 x 10"^ GeV"^ and a^^^,^ ~ 1/137 is the usual electromagnetic fine structure 
constant. In chiral perturbation theory with terms of order up to p^, the amplitudes A and 
B take the form [O]: 
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where we use the kinematic variables 



(91 + ^2) 
Ml 



Pk ■ [qi - q2] 
Ml 



(3) 



and ^ A-nf^ ^ 111 GeV. 



This form for the two amplitudes does not correspond to a complete calculation in 
chiral perturbation theory at order p^. It contains the complete one- loop calculation of order 
and two types of terms of order p^ . The first type consists of the non- analytic terms 
in Eq. |^ that multiply the factors 02 and ai{z). The inclusion of these terms is inspired by 



dispersion relations, and they originate in corrections to the K ^ 3tt amphtudes |n, 
The relevant constants which enter ai and 02 are extracted from an analysis of K ^ data. 
The second type of term consists of the analytic terms that arise from tree-level contributions 
from order chiral Lagrangians. These contributions can be grouped into three unknown 
constants: a^, and j3 corresponding to the three possible Lorentz invariant forms which 
occur at order for the K^TT^'-yy vertex |T^. From the analysis oi K ^ "in in Ref. [0], we 
have 
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In the analysis of Ref. |[TT|], which has become standard, the three unknown con- 
stants were fixed in terms of the contribution they receive from vector-meson exchange, 
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supplemented with a minimal subtraction ansatz: 

a^ = —4a 



a2 = 12ay — 0.65, 



p = -Say -0.13, (5) 



and this form has been used, for example, by KTeV to fit their data with ay = —0.72 ± 
0.05±0.06. In Eq. is no longer independent from 0:1^2; therefore it is clear that this ansatz 
introduces model-dependent correlations between the B amplitude (the one responsible for a 
large CP-conserving Kl —>■ 7r°e~e~), and the A amplitude which dominates the Kl —>■ tt^'J'J 
mode, but which does not contribute significantly to Ki —>■ vr^e^e^. 



In Fig. |l| we reproduce the data from Ref. |16| as can be read from their published 



paper. We superimpose on the data the best fit we obtain in terms of the parameter ay as 
a solid line. Our fit gives ay = —0.95 with a x^/dof = 46/27, which corresponds to 

a-^ = 3.8, 
= -12.0, 

(3 = 7.5. (6) 



Notice that our value for ay is not the same value quoted by Ref. [0 because we 
do not have access to the raw data and hence we have not taken into consideration any 
background or detector issues. Nevertheless, we feel that it is fair to compare this fit to our 
best three-parameter fit obtained in the same way. This one is presented in Fig. |l] as the 
dashed line, and corresponds to 
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Figure 1: Two different fits to the data from Ref. [1^, as explained in tfie text. Tfie solid 
line is a one-parameter fit corresponding to Eq. the dashed line is the three-parameter fit 
shown in Eq. ^ 
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p 



1.7, 
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For this fit we obtain a x^/dof = 37/25, slightly better than Eq. ^. Clearly it is up to the 
experimentalists to present a complete best fit to the data using the general form, Eqs. |I], ^ 

and taking into consideration all the experimental issue^. However, it should be clear from 

^We proceed keeping the branching ratio fixed to the one measured by KTeV (1.68 ± 0.07 ± 0.08) x 
10~^, in the normalization of our fits. Eventually the parameters extracted from all our fits yield branching 
ratios very close to the experimental one and well within its errors. 
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Fig. |I] that even though the current data is consistent with the vector dominance assumption, 
it cannot rule out other scenarios. In fact, our three parameter best fit is not consistent with 
the vector meson dominance assumption. In the Appendix we explore the significance of our 
fit showing the range allowed for its three parameters within one sigma from our best x^- 
This provides an estimate of the errors involved. 

Although the two types of fit are indistinguishable as far as describing the —* 7r°77 
spectrum, they result in completely different predictions for the unitarity bound on the CP- 
conserving contribution to Kl 7r°e"^e~. To evaluate it we need to calculate the absorptive 
contribution from the on-shell two-photon intermediate state to Kl TT^e~^e~, as depicted 
in Fig. 1^. 




Figure 2: Contribution from the on-shell two-photon intermediate state to Bcp{Kl 
7rOe+e-). 



This yields the following bounds on the CP-conserving part of Bcp{Kl ix^e^e ): 



Bcp{Kl ^ 7r°e+e-) > 



2.3 X 10-12 

3.4 X IQ-i^ three-parameter fit 



The above contribution is not the full absorptive part since there is a further cut due 
to on-shell pions. Moreover, the full CP-conserving amplitude includes a contribution from 
the dispersive part of the amplitude, with off-shell photons (and pions). The general form 
of the amplitude is 

Mcp{Kl 7r°e+e") = Gsal^^PK ■ {K+ - ke-){PK + PnT^-ff^v, (9) 

where K is the result of the loop calculation and the extra antisymmetry under fce+ ^ ^e- 
is a reflection of the properties under a CP transformation. Introducing a form factor to 
regularize the virtual photon couplings, an expression for is obtained: 
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(10) 



where s = {ke+ + ke-)^ ■ The log factor is of course expected, since the photon absorptive 
part comes from the expansion log(— s) = logs + m. This representation of the amplitude 
leads to CP-conserving branching ratios: 



BcpiKL 7r°e+e-) = {^-^^ 

I 7.3 X 10~^^ three-parameter fit 



3 Resonance Models for ai, 0.2 and /3 

In this section we present the contributions of scalar and tensor mesons to the parameters 
a^, ^2 1^- ^6 '^ill be able to show that the tensor meson /2(1270), in particular, can 
contribute at a level comparable to that of vector mesons and yet produce a different pattern 
for the three constants. 



We have chosen to follow the notation of |TT|, where the parameters a^, ctg ^^^1 [3 are 



defined by the expression for the invariant amplitudes A and B as in Eq. 0. It is convenient 
to relate these parameters to the three Lorentz invariant couplings that can be derived from 
a chiral Lagrangian at order p^. Writing these couplings as 

C = [c.Klt^'F^^F,, + ^d'^KLd^TT^F^'^F,, + ^d^KLd^n'F^^F^f,^ , (12) 

where F^^^ is the usual electromagnetic field strength tensor, one finds that 

= -2C2 + — , 

C3 

^2 = 4ci + 2C2 + — , 

(3 = -C3. (13) 

The couplings that occur at order in a vector meson dominance model have been 
obtained in [jl2l. They are of the form 



Cv = ^^^{d-K.d^Tr'F'^'^F,^ + 2d^K,d^7.'F-^F,p) (14) 



and, therefore, the prediction of vector meson dominance is that 

= -4ay, 
"2 = 12ay, 

(3 = -Say. (15) 

This prediction is at the heart of Eq. ^, and differs from it only by small additional constants 
which appear in a particular regularization scheme for the loop amplitudes |Tl|. Although 
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this pattern is a firm prediction of vector meson dominance models, a specific value for ay 
is not. For example, in Ref. |T2| the values ay = 0.32 or ay = —0.32 can be obtained 



depending on whether one uses the so called "weak deformation model" or not. This is just 
another way of saying that the concept of "vector meson dominance" is not uniquely defined 
for the weak interactions. In addition, phenomenological treatments of vector mesons such 



as those of Ref. include effects from rj — rj' mixing, which are formally of higher order, 
but which result in significantly different "vector meson" contributions to — > 7r°77. It is 



worth mentioning that a quark model estimate of the parameters a^, ctj and (3 yields 



the same pattern as in Eq. |T5| with ay = {Nc/27)g\{M^/m'^) in the notation of WE 



In addition to the vector meson exchange contributions, the parameters ai, 02 and 
f] may receive contributions from the exchange of scalar and tensor resonances. The effect 
of scalar resonances near 1 GeV has been found to be small [|19|, and we include it here for 
completeness. Moreover, we sidestep the issue of a possible scalar resonance in the vicinity 
of 500 MeV because the physics of this broad enhancement in the J = / = tttt scattering 
amplitude is, to a large extent, already included in the treatment of the pion loops. We 
concentrate instead in resonances near 1 GeV such as the /o(980), and take the simplest 
form for the scalar-pion and scalar-photon interactions (we use U as in the notation of 
Gasser and Leutwyler [pi]]): 



Cs = g.STT[D^UD^U^j + -^g^SF>^- F^,. (16) 

We have not included a coupling of the scalar field proportional to light quark masses be- 
cause it does not contribute to —>■ 71^77, and because there is not enough experimental 
information on scalar-meson decays to extract it. 
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The coupling g^^ can be determined from the decay width of the scalar into two pions. 
Adding the charged and neutral modes we obtain 

nS - nn) = ^Jl-Arl,glMl(l - 2r^, + 4r^,), (17) 

with Ttts = Mj^/Ms- If we identify the scalar meson with the /o(980), and use the particle 
data book figures 5(/o ^ tt+tt-) = 2/3, 5(/o ^ 7r°7r°) = 1/3, [0] and the NOMAD result 
r(/o) = 35± 12 MeV we find g.,^ ~ ±5 MeV (we cannot decide the sign ambiguity from 
the experimental rates). 

The width for the scalar- meson decay into two photons allows us to determine g^. 
We find for the width 



r(s-77) = (^^)- 



47r 



If again we identify the scalar with the /o(980) and use the particle data book value r(/o — > 
7^) = 0.39t[];l|j X 10-3 MeV |2[, we find g^ ~ ±3.9 x lO^^ MeV^^ 

Collecting these results we finally obtain for the contribution of the scalar /o(980) to 
Kl ^ 71^77 (see Fig. |): 

«i = -«2 = -^Gg^g-y^ ~ ±0.08, P = 0. (19) 

In a similar manner we can determine the contribution from a tensor meson. A 
simple look at the low energy data for the reaction 77 — * ti^ti'^ suffices to motivate the 



potential importance of the /2(1270) for our amplitudes through diagrams such as those in 



Fig. H. Following Ref. pO] we write the lowest order couplings of a tensor meson T^^, to pions 
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Kl 



^0 



a) b) 



Figure 3: Scalar- and tensor- meson resonance Feynman diagrams contributing to —>■ 
7r°77. The dots in a) and b) represent flavor- changing mass-insertions in the incoming and 
outcomimg particles, respectively 3; PH- 

and photons as 

£t = KT'^'TtI^D^UDM^^ + ^h^T'^'F^^F^. (20) 



For the inclusive width of the tensor meson into two pions, and following Ref. for the 
description of the spin 2 states, we obtain 

r(T-vrvr) = ^^(l-4m^/M^) . (21) 

For the decay width of the tensor meson into two photons we find 
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Identifying the tensor meson with the /2(1270) and using the particle data book values for 
mass and partial widths p2[, we obtain /i^ ~ ±40 MeV and ~ ±0.03 MeV~^. 



The tensor (/2) contribution to the parameters a^^, can be read from the 

interaction that results after the tensor meson has been integrated out 



47r Mf, V 3 



The resulting contributions are: 



V/3 



(23) 
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Ml 



-Kh,'-^ ±0.25, 



28 ^ Ml 



3 "^""^'M| 



]Vf2 

8KKj^ ~ ±1.5. 



(24) 



We summarize our results in Table |I|. 





Vector {av = ± 0.32) 


Scalar 


Tensor 


Our Best Fit 


Best Fit ay 




T 1.2 


T 0.08 


± 0.25 





3.8 




± 3.6 


± 0.08 


T 1.7 


1.7 


-12 


P 


T 2.4 





± 1.5 


-5 


7.5 



Table 1: A comparison of parameters for —>■ ir^'j'-f for various contributions discussed in 
the text. We contrast these contributions with our best three-parameter fit, as well as with 
our best fit within the VMD ansatz. 



4 Conclusion 



We expect new data for —>■ n^'-f'-f from KTeV and NA48 in the near future, and this 
makes a reanalysis of this mode timely. We have argued that the new results should not be 

13 



analyzed in terms of the vector meson dominance ansatz, but rather in a model independent 
way, and that this entails the use of three parameters. These three parameters are related 
to the three a priori undetermined counterterms entering the amplitude, as shown in Eqs. 

To illustrate the previous point we have re-examined the fit to the 1999 KTeV data. 
We find that the general, three-parameter fit is slightly better than the old fit in terms of 
ay, and we show our results in Fig. |I|. The difference between the two procedures appears 
to be small in the — > 71^77 spectrum. Nevertheless, it leads to significantly different 
predictions for the CP-conserving component of Kl vr^e^e", which can be seen in Eqs. |], 
pJ^ New data, with higher statistics, should be able to better distinguish the two cases. 

As a further motivation for abandoning the usual parameterization, we have also 
shown that the /2(1270) tensor meson can yield an important contribution to the counter- 
terms, and that this contribution cannot be cast in the one-parameter framework of vector 
meson dominance. 
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Appendix 



To compute an estimate of the errors involved in our fits we calculate the range of variation 
of the three parameters from our "best fit" xLin — 37 to xLm + ^ (corresponding to one 
standard deviation), obtaining: 

-2.0< < 1.9, 
0.8 < ^2 < 2.5, 
-5.3 < p < -4.5. (25) 

The "central" values give roughly a linear equation: 

02 + 0.29 = 1.65. (26) 

This is consistent with 2; — ~ 0.3 dominating any z dependence. 
The above values are to be compared with the ones given in Eq. |^: 

ai = 0, 
0^2 = 1.7, 

P = -5. (27) 

In Fig. 1^ we present a one-sigma plot of the parameter space for ai and 0:2 with a fixed 
P = -5.0. 

It is possible to redo the calculation keeping each time one parameter variable and 
the other two fixed at the values of Eq. ^. In this case they are much more constrained: 
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Figure 4: Scatter plot of the parameter space allowed for ai, a-i with a fixed /3 = —5, within 
one sigma from our x^in- 



- 0.5 < < 0.6, 
1.5 < < 1.8, 
-5.1 < /5 < -4.9. (28) 

There exists another region in the parameter space where is within one sigma from 



Xmin" 



- 0.6 < < 0.8, 
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-12.2 < ^2 < -11-8, 



^.4 < f3 < 8.7. 



(29) 



Note that the one-parameter fit in terms of ay lies closer to this second region. 

Fig. ^ is analogous to Fig. ^ assuming f3 = 8.55, and is consistent with the "central 
values" linear equation 

^2 + 0.28^1 = -11.9. (30) 
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Figure 5: Same as Fig. with j3 = 8.55. 
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